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Abstract. We consider the evolution of ultra-short optical pulses in linear 
and nonlinear media. For the linear case, we first show that the initial- 
boundary value problem for Maxwell's equations in which a pulse is injected 
into a quiescent medium at the left endpoint can be approximated by a lin- 
ear wave equation which can then be further reduced to the linear short-pulse 
equation. A rigorous proof is given that the solution of the short pulse equa- 
tion stays close to the solutions of the original wave equation over the time 
scales expected from the multiple scales derivation of the short pulse equation. 
For the nonlinear case we compare the predictions of the traditional nonlinear 
Schrodinger equation (NLSE) approximation which those of the short pulse 
equation (SPE). We show that both equations can be derived from Maxwell's 
equations using the renormalization group method, thus bringing out the con- 
trasting scales. The numerical comparison of both equations to Maxwell's 
equations shows clearly that as the pulse length shortens, the NLSE approxi- 
mation becomes steadily less accurate while the short pulse equation provides 
a better and better approximation. 



1. Introduction 

The standard model for describing propagation of pulses in nonlinear Maxwell's 
equations is the cubic nonlinear Schrodinger equation (NLSE) ^ |2] . Two main 
assumptions are made in the derivation of the NLSE from Maxwell's equations: 
First, it is assumed that the response of the material attains a quasi-steady-state 
and second that the pulse width is large in comparison to the oscillations of the 
carrier frequency 

In most of the applications of the NLSE in the past, i.e., in the case of pulse 
propagation in optical fibers, both assumptions were well satisfied. At present, 
however, technology for creating very short pulses has advanced a lot and experi- 
ments with pulses which are as short as a few cycles of the carrier wave have become 
possible 4, . The description of those pulses lies beyond the slowly varying envelope 
approximation leading to the NLSE [S] . Various approaches have been proposed to 
replace the NLSE - see, for example [SIEIIS] for a sample of these methods. In [5], 
building on work of jlU) . two of us proposed an alternative model to approximate 
the evolution of very short pulses in nonlinear media. 

In the present paper we study further the short pulse equation (SPE) derived in 
[S] . There are two main sections of this paper. In the first we concentrate on giving 
a rigorous justification of several of the assumptions and approximations made in 
PI, for the linearized short pulse equation. Among the approximations used in |H] 
were: 
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(i) The linearized polarizability of the medium could be approximated by the 
Fourier transformed expression 



(see equation (3) of Ref. Typically, in deriving ^ one assumes that the 
medium has reached some quasi-stationary state, and as J. Ranch pointed 
out to us, it is not clear that for these very short pulses the medium will 
have time to reach such a state before the pulse passes, 
(ii) If the expression is an accurate approximation to the polarizability of 
the medium, then can one really approximate solutions of the resulting 
equation which correspond to "short" pulses by solutions of the "short 
pulse" equation derived in [H] and if so, over what time interval does it 
provide an accurate approximation to the true evolution? 
Note that the first of these points is a question just about the linear problem 
and thus in the next section we consider points (i) and (ii) in the context of a 
medium whose polarization is assumed to be linear. We study solutions of an 
initial-boundary value problem in which a linear wave equation is coupled to a 
medium whose polarization is modeled by a damped, linear oscillator. We inject 
a (short) pulse into the left end of this material and study how that pulse evolves 
with time. We prove rigorously that one can, even in this short pulse regime, 
approximate the polization of the material by the quasi-stationary approximation 
([T)l. and we show that if the pulse length is measured by the small parameter e, 
then the linearized version of the short pulse equation accurately describes the true 
solution of the equation over time scales of 0{l/e). 

The second part of the paper studies of the effect of pulse length on the propaga- 
tion of pulses in nonlinear materials. Using the renormalization group method, we 
derive both the NLSE and the SPE from Maxwell's equations. Then we compare 
the (numerical) evolution of solutions of a one-dimensional, nonlinear version of 
Maxwell's equations with approximations given both by the NLSE and the short 
pulse equation. Wc find as expected from the assumptions that underlie the for- 
mal derivation of these two equations that for slowly modulated pulses the NLSE 
does a better job of approximating the evolution, but as the pulses become shorter 
and shorter the NLSE approximation breaks down and the SPE provides a better 
approximation to the true dynamics. 

2. Formulation of the linear initial-boundary value problem 

In this section we study the propagation of a short pulse injected at one end of 
a semiinfinite fiber with linear polarizability. If we assume that the polarization 
of the electric field is transverse to the fiber then the magnitude of the electric 
field, u{x, t), in appropriately non-dimensional units satisfies the partial differential 
equation 



where p{x, t) stands for the (magnitude of the) polarization of the material, which 
we assume is parallel to the electric field. 

We model the polarization of the material by a damped harmonic oscillator, so 
p{x, t) satisfies the equation 




(1) 



Uxx = Utt + Ptt 



(2) 



Ptt + Ppt -I- ujIp = xqu. 



(3) 
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Remark 2.1. In general, the polarization of the medium is modeled as a sum of 
oscillators, each with its own natural frequency and damping constant. However, as 
was remarked in [S], for infrared pulses in silica fibers, the polarizability can be accu- 
rately modeled with a single resonance. Furthermore, in the present circumstances, 
the linear nature of the problem means that if the polarization was modelled as a 
.sum of several oscillators, we could write the solution as a sum of the solutions to 
problems of the type with different resonant frequencies and damping constants. 

Remark 2.2. Experimentally, the damping of these oscillators is quite weak - T 
is small. Thus, throughout this paper we will assume that — AlUq < and also 
that xo > 0. 

Physically, we are interested in the situation where our medium is semi-infinite 
with one end at the origin and all fields in the medium are initially zero. We enforce 
this condition by taking initial conditions 

u{x, 0) = ut{x, 0) = p{x, 0) = pt{x, 0) = for X > (4) 

We inject a pulse into the left end of the fiber and model this by assuming a 
boundary condition at a; = of the form 

u{0,t)^Uo{t) (5) 

where Uolt) represents the optical pulse. Since we are interested in short pulses, we 
will assume that the injected pulse has the form Ua{t) ^Uo{t/e). We also assume 
that Uq smooth and that it is supported in the interval [0, 1]. 

Summing up, the equations |(2Jl, (jSl together with the boundary condition |SJl and 
the initial conditions l@J form the initial-boundary value problem (IBVP) describing 
the short-pulse propagation. 

Remark 2.3. In general one can also expect an instantaneous contribution to the 
polarization - we ignore this case here since it can be incorporated simply by chang- 
ing the coefficient in front of the term Uu . 

We prove two approximation results about the solutions of the IBVP (jSJ . Recall 
that from the formal derivation of the short pulse equation in [H], we expect that 
approximation of solutions of (O by the (linearized) short pulse equation should be 
valid for times 0{-). In fact, following the usual convention in nonlinear optics the 
multiple scale expansion in [5] was made in terms of a long space scale, ea;, rather 
than a long time scale ei, and thus one might expect an approximation result valid 
over length scales of 0(l/e). However, since pulses in Eqn. |(5J) travel with a speed 
0(1), we can translate our approximation result into a result valid over long time 
scales. Our first result shows that for times of this order we can approximate 
solutions of 101 by solutions of the single equation 

Uxx = Uu + Xou , (6) 

with initial conditions u{x,0) = Mt(a;,0) = and boundary condition u{0,t) = 
Uo(t). Note that this corresponds precisely to the approximation in equation (5) in 
[S| and thus gives a rigorous justification of the heuristic argument of that paper. 
More precisely we prove the following: 

Proposition 2.4. Let Tq > be fixed. Then there exist eo > and Cq > such 
that if < e < eo, and u{x,t) satisfies the IBVP 0), |3|), ^} and u{x,t) 
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satisfies with zero initial data and the boundary condition ^ then 

sup (sup \u{x, t)-uix,t)\) <Coe^/^ . (7) 

0<t<To/e x>0 

In this context we can also show that the hnearized version of the short pulse 
equation derived in [H] correctly describes the propagation of solutions in either (O 
or jnj. The linearized short pulse equation describing the evolution of a function 
of two variables U =U{(j), T) is 

2d^dTll = XoU . (8) 

We also prove 

Proposition 2.5. Let Tq > be fixed. Then there exist eg > and Cq > such 
that if < e < eq, and u{x,t) satisfies the IBVP {0), 0^, ^} then there 

exists a solution of the pulse equation 0), U{(j),T) such that 

sup {sup\u{x,t)-U{*—^,et)\)<Coe^/^. (9) 

l<t<To/e x>Q e 

Remark 2.6. Note that in this proposition we do not begin to compare the solution 
of the "true" evolution 0) with the solution of the pulse equation until a time t > 1 
- i. e. until the pulse has been injected at the left boundary of the domain. This is 
because we don't expect the pulse equation to describe the evolution before there is 
a pulse present in the system! 

Remark 2.7. In both of these propositions the constants Cq depend on the profile 
of the injected pulse, Uq, in a way we make precise in the proof. 

We note that in 11 , Alterman and Rauch study in detail the properties of a 
linear short pulse equation similar to ijSJ, but lacking the term xo^- 

3. Proofs of the approximation results 

We begin with the proof of Proposition 12.41 In the proofs of both propositions 
we will work with the Fourier sine and cosine transforms. Given a function v{x) 
defined on the positive half-line we define: 



oo 



ii^k) — / sm{kx)v{x)dx , (10) 







oo 



v'^{k) ~ / cos{kx)v{x)dx . (11) 







We will also use the Laplace transform which we denote by £. If we take both 
the Fourier-sine and Laplace transforms of (j^J and use the boundary and initial 
conditions we find 

-k^£[u'] + kC[Uo] = s^C[ii'] + s^C[p'] , (12) 
{s^ + sr + Luf,)C[f]^XoC[u'] . (13) 
We can combine these two expressions to obtain a single equation for namely 



C[u^]ik, s) = { ^ \ C[UoKs) . (14) 
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If we now take the inverse Laplace transform of this expression and use the fact 
that C/o(t) = we find that 

u%k,t) = e f' J -. ^ i - a])U,{a)da . (15) 

We now rewrite this expression with the aid of the following standard lemma about 
Laplace transforms: 

Lemma 3.1. If F{s) =£[/](s), then 



e 



f{et)=C-'[F{-)]{t) 



Using Lemma l3 . 1 1 and defining p = ek, we find that 

*_ 

u^^,t)^ [' £-^[F{s,p;e)]{-~a)Uo{a)d<j , (16) 
£ Jo e 

where 

F(s p- e) = (6p)(,s^+e.r + eH^) 

One can compute the inverse Laplace transform of F via it's partial fraction expan- 
sion and for that we need to find the roots of the polynomial in the denominator 
of F, i.e. 

Q{s;p,e) = (s2 + esT + e^ojl)(p^ + s^) + e\os' • (18) 
To this end we use the following series of Lemmas whose proofs are elementary but 
somewhat involved. Hence we relegate the proofs to Appendix IXI 

Lemma 3.2. For all values of p and e all the eigenvalues of Q have non-positive 
real part. 

Lemma 3.3. There exist eo > and Co, Ci > such that for |e| < eo andp > C^e, 
Q has a pair of roots of the form 



sl = --{T±^T^~Au:l) + eal{p) 

with 

< ^ . 

Lemma 3.4. There exists eg > and Co > such that for |e| < eo o,nd p > Cge, 
Q has a pair of roots of the form 



S± = ±? vV + ?XO + Cr± {p) 

with 



p2 



Furthermore, the real part of a\. is negative. 

Lemma 3.5. Let eo and Co be as in Lemma \'d.4\ There exists Cmm > such that 
forp < C'oe the roots of Q can be written as s±(p) = es\{p/t) and s)^{p) = es)^{ple) 
with s\ (q) and s}j_ {q) all distinct. Furthermore, 

min(|S° (g) - S°_{q)\, \sl{q) - si{q)\) > c™„ , 
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while 

Remark 3.6. Note that a corollary of the vroof of Lemma \S. ,^ is that for p > the 
roots of Q are all distinct. 

With these estimates on the eigenvalues of Q in hand we now construct the partial 
fraction expansion of F. Note that since the eigenvalues of Q depend continuously 
on p (and e) we can label the eigenvalues as s\{p) and sj_(p) for all p > 0, and then 
one can write 

Al{p) A" [p) A\{p) (p) , ^ 

s ~ s^^[p) s — s_[p) s — s^(p) s — s_(p) 

If we also write 

F(s V e) = iep){s^+esT + e^xo) ^20) 

^ ' is - - - s\ip)){s - slip)) ' ^ ' 

then we see that A^ has the form 

+ iis'i)-s°_)iisO^)-sl)iisl)-sl) ■ ^''> 

First note that Hp < Coe, the numerator of this expression can be bounded by Ce^ 
by Lemma l3 . 51 while the same lemma guarantees that the denominator is bounded 
below by ce'^ for some c > 0. Thus for p in this range |A!|_| < Ce. 
Now suppose that p > Coe- Since s!}_ is a root of Q, we have 

\iisl)'+e{sl)r + eW,^' ^'^"(^+)' 



Now we see that this expression is bounded by Ce^/p^. Using the asymptotic 
expressions for the roots of Q coming from Lemma 13.31 and l3 . 41 the denominator in 
l|21|l can be bounded from below by tp^ in this range and hence we can bound 

\K\ < % , 

for p in this range. We can combine the estimates on \AP^\ in these two ranges of p 
along with identical estimates on A^_ to obtain 



Lemma 3.7. There exists a constant C'a > such that for all p > 0, 



p- 

We now estimate the coefficients A}^. Following an argument similar to that 
above we find 

,1 _ (6p)((4)^ + e(4)r + e^a;g) 
+ iis\)-sl)iis\)-sl){{s\)-sO_) ■ ^''> 

For p < Coe, using the asymptotic values of the roots and given in Lemma 
13. 51 we see that 

\Alip)\<Ce. 
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For p> Coe we use Lemma [3.31 and rOI to rewrite 

{sir + + e^u^l {{s\f + e{s\)V + e^u^l) 



where 



Thus, 



((4) - - ^'^-) ((4)' + - « + + 

= l+f(p;e), (23) 



\E{p-e)\<-^^. (24) 



^1 ^ ep(l + g(p;£)) 



or if we write 



A\ = ^ + AA^ , (26) 



we see that for p > Cqe, (and Cq sufhciently large) 

Similar estimates hold for A\_ and we have established: 
Lemma 3.8. There exists a constant Ca > such that for p < Coe, 

\Ai{p)\ < CAe , 

while for p > Cqc, 

A^^ip)^-^^L= + AA^^. 
2iy/p-' + e-^xo 

with 

CAe'p 



\^Aiip)\ < 



(p2 + e2)3/2 ■ 



With these estimates on the coefficients in the partial fraction decomposition of 
F{s,p]e) in hand, we now return to the task of computing the solution u{x,t) of 
(0). From and we see that 

r\ poo /> - 

u{x,t) - — / sin(p(-)) / ' Al{p)e''>^^i-''^Uo{a)d<jdp 



sin(p(-)) r A"_{p)e'-'^P^'^-'-'''>Uo{<j)d(jdp 



TTc Jo e 

(p(-)) / ' Alip)e'+^P'>^i-^^Uoia)dadp (28) 



2 . , ,x 

H / sm 

+ — / sin(p(-)) I Al{p)e'-^P'>^i-'''>Uo{(j)dadp . 
Tre Jo e Jo 

We can immediately bound the first two terms on the right hand side of H28|l by 
using the facts that the real parts of Sj_ are negative, so that the exponential factor 
is bounded by 1, as is the factor of sin(pa:/e), and bounding Aj_ by the bound in 
Lemma 1X71 Integrating over p and a we then see that these two lines are bounded 
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by CellZ^ollii- We now turn to the last two lines in l|28|l. First of all, rewriting them 
with the aid of Lemma [3. 41 and as 



2 

7re 



:sm(, , , 

• e 



in(p(-))sin(\/p^ +e^Xo(- - cj))Uo{a)dadp 



2 

Tre 
2 

7re 



p=CiV7"'0 



+ AAL (p)e^- ( ^ ) sin(p( - ) )Z^o i(T)dadp 

e 



1 



sin(p(-)){ 



ep 



2* + e^Xo 



ep 



(29) 



2« Vp^^^^o 

+ — r ^'"^ /' sin(p(-))(Ai (p)e^+(P)(^-") + A^p)e*-(p)(^-")|z^o(^T)rf'Tdp 

The last of these integrals can be immediately bounded by CA\/e||Z^o||Li ■ The next 
to last integral is estimated by using Lemma and the fact that for < f < To/e, 

Ce 



|2 ■ 



p 

With this estimate the integral can be bounded by 



C 



P 



e'p 



(p2 + e2-)3/2 I p2 



Finally, in the second term we use the bounds on A^^^ from Lemma l3.8l to estimate 
this integral by 



P 



;dp\\Uo\\L^ <Ce\\Uo\\L^ ■ 



Combining these estimates with those on the first two integrals in H28|) we see that 

2 



sup sup 

0<t<To/e \ 2;>0 



u{x, t) 



'\ sin(p(-))x (30) 



X sin ( v'p^ + e^Xo(- " cr) ] Uo{a)dadp 



We now note that since 

2 /-^^v^ n ep 



: sm 



(p(-))sin(v/p^ + e^Xo(- - (T))Uo{a)dadp 



'p=0 Jo ^p2 + e2^jj 

< C^/i|j^^o||Ll (31) 

we can subtract it from the left hand side of 1)30(1 without changing the bound on 
the right hand side of this expression. That is, we can bound the difference between 
u{x, t) and the integral 

2 / " / » ep 

^=s^np^-jjs^nvP" + e-Xol^ 



P=o Jo ^p2 + e'^xo 



X t 

(p(-)) sin( Vp2 + e2xo(- - a))Uo{a)dadp. 
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by C^UoWl. 

By taking the Fourier sine transform of ijHJl we see that this integral is exactly 
u{x,t), the solution of © and thus we have completed the proof of Proposition 
O □ 

We next prove Proposition 12.51 By the results established so far it suffices to 
prove that the solution t2(x, t) of © can be approximated by a solution of the pulse 
equation ^ over the relevant time intervals. First note that using trigonometric 
identities for the sine and cosine we can rewrite 



u{x,t) =u^{x,t) +u'^{x,t) (32) 



where 



u^(x,t) ^ — / / — — cos(-( px - uj J p)t) + auJe{p)Wo{o')dadp , (33) 

TreJo Jo ^Ap) e 

and 

u^(x,t) = / / — —coii{-{px + uj^{p)t)-auj^{p))Uf){a)dadp, (34) 

where uje{p) — \/p'^ + e^Xo- 

Roughly speaking, and represent the left and right moving parts of the 
pulse. In particular, for t > 1, we expect that the left moving part of the solution 
will no longer be relevant since we are only interested in the solution for x > 0. To 
prove this we use another trigonometric identity to rewrite as 

1 p 1 

u^{x,t) = / — —cos{-{px + Ue{p)t)UQ{uj^{p))dp 

TT Jo We(p) e 

1 p 1 

■sm{-{px + uJe{p)t)UQ{uJe{p))dp , (35) 



TT Jo ^eip) e 

and we recall that Uq and Uq are the cosine and sine transforms of the boundary 
data Uo- (We have used here the fact that since the limit on the cr integral exceeds 
the limits on the support of Uo we can integrate from to oo.) We now prove that 
both of these terms are C(\/e) in the L°° norm and thus can be ignored to the 
order of approximation that we are concerned with. 

We'll consider in detail the first of the two terms in H35|) . The second is handled 
in an almost identical fashion and we leave the details as an exercise. Rewrite that 
integral using a trigonometric identity as 

1 p t P - 

/ {cos{-ix + t))cos{-{uj^{p) ~ p)}——U^{uj,(p))dp 

n Jq e e ujAp) 

I {smi^ix + t))sm{-iLu,ip)-p)}-^^U^{co,ip))dp . (36) 



Once again these two integrals are estimated in an almost identical fashion so we 
provide the details for the first and leave the second as an exercise. Integrating by 
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parts, the first integral becomes 



1 e 

■ sini 



{l(:, + t)){iU{p) - l)-^sin(-(^e(p) -p)M(c.e(p)) 



n Jq X + 1 e ^ e ^ (p) 

-cos(^(.;,(p)-p))[ ^-^^^~^^y^^ K(c..(p) (37) 

-cos(-(a;,(p) -p))^a;^(pK'(c.,(p))}dp . 

Note that there exists a constant Ci > 0, independent of e such that the the various 
quotients appearing in the integrand of H37I) can be bounded as follows: 



jCi forallp>0, 



^e{p) ' " forallp>Co^/i. 



, cj,(p)~pa;^(p) , ^ j forallp>0, 
' '"l^ forallp>CoV^. 



iP'^'eip) I 



< Ci for all p > 



(|i| + l)K(c.,(p))|dp (38) 



Thus, bounding the factors of sine and cosine in H37|l by 1 we see that this integral 
can be bounded by 

x + t Jo 

+—7 / (1^1 + 1)(- + -)\mi^.mdp+—-. / \w,'{i^M)\dp- 

In the first two of these integrals we bound \U§{Ljf{p))\ by C||Z^o||li Thus, these 
first two integrals can be bounded by Cy/e\\Uo\\Li. To bound the final integral write 
it as the sum 

f 1 poo 

\Wo'{^e{p))\dp+ / \U^'{uj,{p))\dp 



The first integral can again be bounded by C||Z//o||li, while the second is bounded 
by C /|°° |Wo'(^)|d^, by making the change of variables ^ = We(p). This integral can 
be bounded by C{J^{1 + ^^)\U§'{0\^dO^^^ by the Cauchy-Schwartz. Applying 
Parseval's equality and the fact that Uq has finite support this integral is bounded 
by ||i^olii/i ■ Note that since Z^o has finite support, one can also bound the norm of 
1^0 by a constant times the norm, and thus, Ip^ is bounded by ^/e\\Uo\\H^ 
A similar estimate applies to the remaining terms in and so we conclude that 

snpsup\ii~u"'{x,t)\<C^/7\\Un\\m ■ (39) 

t>l a:>0 

We now examine more closely and show that it can be approximated by a 
solution of the pulse equation. Begin by writing it as 



u ^{x, t) ^ - — --- cos{-{px - uj^{p)t)UQ{uj^{p))dp 

' ^ ■ sm{- {px ~ Lu,{p)t)U^iLj,{p))dp . (40) 
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Now define (p = {^) and T = et. Then 

z^«(0,r) = 7.«(x(0,r),i(0,r)) 

1 p T 



Jo ^e{p) 



cos{p(j3 +{p- uj,{p))^)U^{u:,{p))dp (41) 



We now define 



1 ^ 
+ - sm{p(j> +{p-u;, {p) ) - (we {p))dp 



TT J ^1/2 U)e{P) 2p 



P .sin(p0-^K^(c.e(p))dp. (42) 



TT LUeip) 2p 

Note that by an easy and exphcit computation W(0, T) satisfies the linearized pulse 
equation (jSl, hence, Proposition 12 . 51 will follow if we can show that AU ~ — U 
is small. Subtracting H42|) from (|41|) we obtain 

AZY(</., T)^- [ |cos(p</. + (p - ^.(p))^) - cos(p0 - ^)\w^{u;,{p))dp 

H — / — -- <^ sm{p4) +ip- uje{p))^) - sin(p0 - ^^) I U^{uj,{p))dp 

H / — ^ COs(p(/) + (p - U}e{p))^)ilo{uJe{p))dp 

Jo ^e[P) 

+ - / ^ sin(p0 + (p - u;,{p))^)U^{u;,{p))dp . (43) 

Note that the integrals over [0, y/e] can be easily bounded by noting that jp/cj, (p)| < 
1, hence both of these integrals are bounded by C^/e(|Z^Q + The two 

remaining integrals are bounded in a similar fashion - we give the details of the 
bound on the integral containing the difference of cosines and leave the other as an 
exercise. By the mean value theorem there exists some ^ such that 

\ cos{p(f> + {p - uj,{p))^) - cos{p(f> - ^^)| = |sin(^) ((p- ^^(p))^ + 
<CT^ , 

p6 



— j| = |sm^^j^^p-a;HPJJ^ + ^J 



where the last inequality used Taylor's theorem to bound | (p — (p) ) ^ + 

and the fact that for p > e^^^, ^ << 1. Inserting this estimate into the first of the 
integral terms in (|43|l we see that it is bounded by 

7,1/2 p'^UJeip) 

A similar estimate holds for the remaining term in the definition of AU and Propo- 
sition |^| follows. □ 
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4. Approximating the nonlinear pulse dynamics using renormalization 

GROUPS 

Summarizing the results of the previous section we now know (rigorously) that if 
we ignore nonlinear effects we can approximate the motion of a short pulse injected 
into one end of an optical fiber by the linearized short- pulse equation 0. Since 
nonlinear effects are important for many optical phenomena, the next natural step 
is to investigate how incorporating nonlinear terms into the polarization affects (jSJ . 
Here, we consider as a first step the simplest form of a nonlinear contribution p^i 
given by 

Pn\ = X3U^- (44) 
In order to answer this question we start from a wave equation similar to J^Jl with 
an additional nonlinear term 

Uxx ^ utt + Xqu + X3iu^)tt- (45) 

Remark 4.1. The equation \45}) can be derived from Maxwell's wave equation 

Uxx - Utt = (Plin)tt + (Pnl)tt 

writing the linear part of the polarization as 



Plin{x,t)^ J X{t - T)u{x,T)dT 

and making the approximation 

This means physically that the frequency range of the pulse under consideration 
is not only far from the resonance frequency of the material but also much larger 
than loq. It is also possible to consider other forms of the susceptibility in frequency 
domain, leading to different types of equations for the short pulse. If we ignore the 
nonlinear term, then Proposition \2.4\ shows that this approximation leads to a small 
error. The choice of nonlinear part of the susceptibility corresponds to 

Pni{x, = J x'^^ {t - Ti,t ~ T2,t ~ T3)u{x, Ti)u{x, T2)u{x, Ta ) ^Ti (ir2 dra 

with the assumption that the nonlinear contribution is instantaneous, hence 

X^^Ht - Ti,t - T2,t- Tg) = X3S{t - n)5{t - T2)S{t - T3). 

On the basis of formal asymptotic calculations we believe that this is the most im- 
portant contribution to the nonlinearity and thus in the present paper we limit our- 
selves to the consideration of this case. However, we also stress that, for ultra-short 
pulses, it is very interesting to extend the analysis to more complicated forms of the 
nonlinear susceptibility including delay in the response of the material. 

As mentioned in the introduction, the standard model describing the nonlinear 
pulse evolution is the cubic nonlinear Schrodinger equation (NLSE). To emphasize 
the different regimes to which the NLSE and short pulse equations apply, we briefly 
review how one derives the NLSE from (|45|1 . The main idea is that we assume a 
broad, rather than a short, pulse in the sense that we introduce time scales that 
are slower than the oscillations of the carrier wave that is oscillating at a fixed 
frequency Cj with a wavenumber (3. Therefore, the NLSE is an equation describing 
the slowly varying amplitude of the optical signal. 
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The separation of those time scales can be done by a usual expansion in multiple 
scales 0. In the present work, however, we utilize the so-called renormalization 
group (RG) method to derive the NLSE. This perturbative technique was first 
developed by Chen, Goldenfeld, and Oono as a tool for asymptotic analysis (see 
and [n]). In [HI, the validity of the RG method has been justified by applying 
it to various examples of ordinary differential equations involving multiple scales, 
boundary layers and WKB analysis. See also, jTH], JEli and jTZj for some examples 
of the rigorous use of the renormalization group in the study of partial differential 
equations. The mathematical study of this method has also been presented by 
Ziane in ^H] . The author explicitly described the RG method in the general setting 
of autonomous nonlinear systems of differential equations. 

We will follow the approach given in I18J and explain how to obtain the per- 
turbative solution of Maxwell's equation. Although the RG method will lead to 
the same results as in multiple scale technique, it is worth mentioning that there 
are advantages of using this method. First, the RG method does not require one 
to introduce all the different scales in the beginning of the ansatz since these will 
appear naturally in the renormalization group equation. This implies that one can 
assume a naive perturbation series in any given problems involving multiple scales. 
A second argument in favor of this method is that the algebraic calculations are 
simpler than when other perturbation techniques are used, especially when one 
considers higher order approximations. 

To see how the RG method works in the present case, we start from a slightly 
more general form than (|45|l 

Uxx = uu + dtt J x{t - T)u{x,'T)dT + X3{u^)tt- (46) 

We assume that the solution of (|46|l is of small amplitude and concentrated around 
the carrier frequency. Because of the oscillations of the carrier wave, in the Fourier 
domain the signal will be concentrated around the frequencies Co and —uj. Therefore, 
we can write our solution as a wave packet in the form of 

u{x,t) = C/(a;,i)e*"^^~"*) +C/*(a;,i)e"'^^^"^*^ (47) 
Taking the Fourier transform of H46I) . we find 

+ ^(^' ^) = -^^X3^{x, iu), (48) 

where the wavenumber /3 is given by 

(3{uj) = uj^l + x{^). (49) 

The main idea is now to make a Taylor expansion of the dispersion /? around the 
carrier frequency. This assumes that the signal is localized in Fourier domain, cor- 
responding to a slowly varying amplitude approximation in time domain. Because 
of this local character of this expansion, the specific form of xi'-^) is not essential. 
The Taylor expansion of /? at (2; yields 

{Lu-Hjf + ■■■ . (50) 



80' 
duj 



{lj — Oj) + 



2 diu^ 
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Let US denote /3 — /3(a)). Applying the inverse Fourier transform to iPHji . the 
straightforward calculation yields 

k 



u{x,t) + e'i0^'^t)J2 



1 9*^/32 



.d_ 

dt 



U{x,t) 



^ kl dbj^ 

k 



52 



(51) 



Since we assume that the signal is concentrated around a carrier wave , we now 
introduce a slow time ti by setting 



Then we can write (|47ll as 

u{xMM) = t/(x,ti)e'(^^-'^*«) + C/*(a;,ii)e-*('^^-"*«\ 



(52) 



(53) 



where U{x,ti) is a slowly varying function of x. Here, notice that we do not 
explicitly separate the scales in the evolution variable x at this step. Now we use a 
small amplitude expansion of the function u 



euo + e^ui + e^U2 + ■■■ . (54) 

Solving now 1)51(1 order by order, we can determine the equation for U{x,ti) which 
describes the slowly varying amplitude of the electric field u. 

Without any loss of generality we assume X3 = 1. We first collect terms 0{e) 
and find 

■^+(3^^UQ{x,to,ti)^Q. 

Assuming the solution has the form of 1(53(1 . we plug it into the above equation. 
Due to the ansatz ((53(1 . at the lowest order we find the amplitude of solution uq 
does not depend on x and obtain 



Mo 



(x,io,<i) = Ao(ti)e*('^"-'^*") +^^(ti)e~^('^""'^*°). 



Here, ^0(^1) can be explicitly determined from the initial condition for ((46(1 . In 
fact, one can verify for the homogeneous equation (-r-^ + p)u ~ 0, the amplitude 

OX'' 

of the solution does not depend on x due to our ansatz. We will now see how the 
a;— dependence of amplitude enters. At the second order, ©(e^), we find 



(5^+/3^j-i(-^*o,^i) + ^ 

+ ^ 



i-^Al(tx)\ e-»(/9--"*o) ^ 0. (55) 



We will write the solution of 1(55(1 as the sum of a particular solution P{x,to,ti) 
and the general solution ^1 (ii)e*('^^-'^*o) + (ii)e-*(^^-'^*o) of the homogeneous 
equation, i.e., we write 

ui{x, to, ii) = ^i(ti)e^(''"-^*") + (ti)e-^(^"-^*") + P{x, to, h), 
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where Ai{ti) depends on the initial condition and we denote P{x,to,ti) by the 
particular solution to equation H55|l . One of the simple ways to find this particular 
solution for the given equation is by letting 

P{x, to,ti) = e'^^'aitQ, ti)x + e-'^'^bito, ti)x, 

where a{to,ti), b{to, ti) are to be determined later. Plugging this into equation l|55|> 
yields 

\ dti dti 
Then, we obtain the second order approximate solution 
u''^\x,ta,ti) = e(Mo + eui) 

+ complex conjugate. 
Letting Ao{ti) = Ao{ti) + eAi{ti), we find 



u 



+ complex conjugate. 

Since u^'^\x, to, ti) needs to be an approximation valid up to order ©(e^), the term 
0{e^) can be neglected. Hence, we have 

u^'HxM) - e(^Ao{h)~eP'^Mll,y^(A^^^to) 

+ complex conjugate. 

Here, we notice that a secular term appears, which corresponds to the term pro- 
portional to X. In other words, this approximation is no longer valid when x 
O (1/e) or higher. In order to get rid of this secular term, we consider the term 

Ao{ti) — ef3' — y^~^x as the order 1 Taylor expansion of some function A{x,ti) 
dti 

about a; = 0. Thus, we need to find A(x,ti) which satisfies that 

A{x,h)\^=o = (56a) 
dA{x,h) ~^, dA{xM) 

This is renormalization group equation. Now the above form of the equation mo- 
tivates us to introduce a new scale ex. Let us define xi = ex then H56|l gives 

A(a;i,ii)U,=o = Mh), (57a) 
dA{xi,ti) ^ ~, dA{xi,ti) 
dxi dti 
By solving Ij57bp provided that the initial condition 157al) is satisfied, we can express 

Ao{ti) — /?'^— 2^— as the Taylor expansion of A(a;i, ti) about xi = 0. Then, we 
dti 
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finally obtain the second order approximate solution without a secular term 

u^'^'> = eA(xi, ii)e'('^'="'^'o^ + complex conjugate. (58) 

In fact, we can find the explicit form of A(xi,ii), that is A{xi,ti) — Ao{l3'xi —ti). 
We notice that the solution of the RG equation A(a;i,ii) remains 0(1). Thus, the 
second order approximate solution is valid over the time scales of interest. In the 
same way we can find the contributions of higher orders. A detailed calculation can 
be found in the Appendix. Since we want to analyze the effect of nonlinearity, we 
need to explore the order 0{e^) where the effects of nonlinearity first occur. The 
result of the application of the RG method yields 

u{x, to, h) = eY,{x, ii)e*(^^~'^*'') + e±*{x, ti)e-'('^^-"*«) + 0{e^) (59) 

where the function S(x,ti) satisfies the NLSE of the form 

df:{x,t,) _ ^~^, dt{xM) (go) 



dx dti 

(-f ^S(a;,ti) + ^is(x,ii)|S(a;,ii)p 



In order to carry out a direct comparison between (|45|l and (|60|) the only remaining 
step is to specify in the dispersion /3(w). We calculate this from the assumed 
form of the susceptability, x{uj) = — xo/^^- In this case, (|49|) yields 

P' - r4 ' ^" = ( 2'''°W2 - (61) 

As mentioned before, we expect the cubic nonlinear Schrodinger equation to be a 
good approximation for fairly broad pulses. For short pulses the scaling is entirely 
different and, based on the results of the previous section based on the linear case 
(PJ, it is reasonable to assume that for short pulses we can approximate the solution 
of (|45|l by an ansatz with the scaling of the form 

u{x,t) = eVQ{^--^,ex) + e^Vi{^—^,ex) + ... (62) 
e e 

into ()45|l . The leading order, nontrivial behavior of Vq is then given by [HI 

-2d^dxVo = XoVo + X^d^^i^of, X = ex. (63) 

This nonlinear short-pulse equation describes the influence of the nonlinear con- 
tribution of the polarization to the pulse. Note that there is a slight difference in 
the nonlinear ansatz (|62|) from the linear result presented in lO, namely that the 
evolution variable is ex in contrast to et. Since our system is weakly dispersive, the 
pulses propagate on the leading order with speed one, an approximation valid to 
0{l/e) in time should also be valid to 0{l/e) in space. 

Eq. (|63|l can be derived by the renormalization group (RG) method as well. 
Starting from H45|) 

Uxx = utt + Xou + X3{u^)tt- 
we first make a coordinate transform 

u{x,t)^B(^—^,x] (64) 
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which takes care of (a) the introduction of the time-scale for ultra-short pulses, 
namely that the initial condition will depend on t/ e and (b) a moving frame corre- 
sponding to a right-moving pulse. The transformed equation for B{(j),x) is 

-^S^x -XoS-B,, + X3 (5^)00 (65) 

As we will see, the attempt to solve this equation by a standard, weakly nonlinear 
perturbation expansion of the form 

B(0, x) ^ eAo{(l), x) + e2Ai(0, x) + ... (66) 

will lead to secular growth: At the leading order, we obtain that Aq does not depend 
on X. This reflects the fact that, in our scaling, the leading order solution is the 
initial data traveling to the right with the speed one. Or, in other words, the system 
we are considering is weakly dispersive and nonlinear such that neither dispersion 
nor nonlinearity enter the solution on the lowest order. At the next order, however, 
both effects are present, since we find as a solution for Ai 

Ai^-ix r xoAo(0') + X3(Ao(0')')00#'- (67) 

^ J -oo 

The expansion for B{(j)^x) is then 

B{<P, x) « e I^Ao - I X XoAo(0') + X3(Ao(<^')')00 ^0'^ ■ (68) 

In order to eliminate the growth of Ai with respect to x, we look for a transformation 
of the form 

Vo(</),x = 0) = Ao(0), (69) 
dVo _ € 
dx ~ 2 , 

or, introducing a slow variable X = ex and setting Vo = Vo(0, A'), we find 

-2dxd^Vo - XoVo + x^dl {VI) (71) 

which is exactly the short-pulse equation H63|) . In this derivation we see again, how 
the scales of the evolution variable x appear naturally. Of course, we could have 
also used a multiple-scale expansion of the form 

B{4),x) = €Va{(i>,xa,xi) + e^Vi{(l),xa,xi) + ... 

which leads, as the reader can easily verify, to the same result. 

We expect on the basis of the derivation of both approximate equations and 
extensive numerical and experimental evidence that for broad pulses the nonlinear 
Schrodinger equation is an excellent approximation but for ultra-short pulses, the 
nonlinear short-pulse equation should be a more appropriate approximation than 
NLSE. Intuitively this is clear by the scaling of the ansatz that was used in both 
derivations: In order to derive (|63|l we started from a pulse of the form U{t/e) 
whereas the NLSE describes the envelope on a slow time scale A{x,et). On the 
other hand, it is not clear how far we can push each of those perturbations. It would 
be interesting to compare (|6;-{|l and H6()|l analytically to the solution of Maxwell's 
equations given by H45|l . but this is an extremely complicated problem. Therefore, 
we approach this problem numerically. 



XoM^') + X3{M4>r)4>4' deb' (70) 
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5. Numerical comparison of the approximations to Maxwell's 

equations 

We perform the following (numerical) experiment: Consider Maxwell's equation 
(|45|l with an initial data that corresponds to pulse carried by a carrier wave: 

u{x = 0, i) = a e"*"'*'/^ cosujot. (72) 

Here, we choose x to be our evolution variable and choose Ux{x = Q,t) such that 
our initial conditions correspond to a forward traveling wave of the linear problem. 
The factor a in (|72(l corresponds to the amplitude of the pulse and the parameter 
b determines the pulse width. For the susceptibility x{^) use in the numerical 
simulations the following form 

x{ij) ^^^iniij-iJc) + ni-iJc-Lj)) (73) 

with 7i (a;) = 1 for x < and 7i (a;) = elsewhere. This accounts for the fact that the 
pulse cannot propagate for very low frequencies. We assume that the amplitude of 
the pulse decays to almost zero at lUc- This condition is satisfied by an appropriate 
choice of the carrier frequency loq > lOc- First, we discuss the "classical" case for 
small 6, where the NLSE applies. To apply the NLSE model, we first extract out of 
(|72|l the corresponding initial condition for the slowly varying amplitude, compute 
the evolution of this initial data according to 1()U|) and then use (|59|) in order to 
reconstruct the electric field u(xcnd,i)- Figure ^ shows the typical result of the 
case where the nonlinear Schrodinger approximation holds. Due to the choice of 
h — 0.2, the width of the pulse is sufficiently large in comparison to the period of the 
carrier frequency. Since we are interested in the nonlinear evolution of the signal, 
we present the result of the NLSE approximation in the following form: First, we 
propagate the initial pulse in the linear Maxwell's equations by setting ^3 = 0- 
This solution itiin serves as reference data: Now we propagate the same pulse in 
the corresponding nonlinear setting and, after obtaining Wmaxwcib we compute the 
difference from the linear solution 

^^maxwell ~ ^maxwell ^lin- ('^4) 

Then we find by solving the NLSE the corresponding Unisc f^nd compute 

AUnlsc = Unlsc ^ U\in- (75) 

The question is now how well Aitniso approximates Aumaxweii- As we can see from 
figure n in this case, the approximation is excellent. Let us now increase b corre- 
sponding to making the pulse shorter. Let's first choose 6 = 2. Notice that this 
parameter already formally violates the basic assumption about separation of time 
scales made in the derivation of the NLSE. Therefore, from figure|21 it is surprising 
how well the NLSE still works. In this parameter regime, the short-pulse equation 
is not better than the NLSE since for 5 = 2 the essential assumption about a t/e de- 
pendence of the initial condition is not satisfied. It is possible to extend the validity 
of NLSE to shorter pulses by incorporating higher order terms. In the Appendix, 
we also give a derivation for the next order that appears in the RG expansion. In 
this article, however, we want to focus on the comparison between the leading order 
approximation of Maxwell's equations and the ultra-short pulse equation. Going 
to shorter pulses, e.g. for h = 3.0, already, the short-pulse equation starts to do a 
better job than the NLSE as we can see from Figure El Setting 6 = 5, we finally 
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Figure 1. Comparison of the solution of Maxwell's equation and the 
cubic nonlinear Schrodinger equation. The figure above shows the initial 
pulse at a; = 0. The figure below compares the difference of nonlinear 
Ma^xwell's equation and the solution of the corresponding linear problem 
to the prediction of the cubic nonlinear Schrodinger equation. The total 
propagation distance is x^nd = 50. The parameters of the simulation 
are a = 0.1, b — 0.2, uj^ = 2.5, uiq = 4, xo = 5 and xa = 0.5. Eq. (|MJ 
yields ^' ?a 1.2. 



arrive to the domain of ultra-short pulses. Here, the NLSE still predicts the rough 
shape of the pulse, but does not give correct information about the pulse shape as 
we can see from figure |21 On the other hand we can see now that the short-pulse 
equation already in this chosen case of 6 ^ 5.0 corresponding to e = 0.2 provides 
an excellent approximation of Maxwell's equation. Note that we propagated till 
a^cnd = 50 '--^ O (^) . This numerical experiment cannot substitute for a more thor- 
ough analytical investigation, but it is an encouraging sign that, for ultra-short 
pulses, can be used in order to approximate the solutions of (|45|l . In further 
numerical studies we found that even at propagation distances up to Xond — 200, 
the short-pulse equation is a good approximation to Maxwell's equations. 
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Figure 2. Comparison of the solution of Maxwell's equation and the 
cubic nonlinear Schrodinger equation and the short-pulse equation for 
a short pulse. Again, the figures compare the difference of nonlinear 
Majcwell's equation and the solution of the corresponding linear problem 
to the prediction of the cubic nonlinear Schrodinger equation and the 
short-pulse equation. The parameters for this simulation are the same 
as in figureQwith the exception of b = 2, uiq — 6.5 for the figure above 
and b = 3, LJo = 13 for the figure below. 



6. Conclusion 

In this paper we presented two main results: First we showed that if we ignore 
nonlinear effects one can rigorously approximate the evolution of a very short pulse 
injected into one end of an optical fiber by a solution of the (linear) short-pulse 
equation (jHJi. Second, in the case of nonlinear pulse propagation, we have derived 
both the cubic nonlinear Schrodinger equation and the short-pulse equation directly 
from Maxwell's equations using renormalization groups. We have shown numeri- 
cally that the (nonlinear) short-pulse equation gives an excellent approximation to 
the solution of Maxwell's equation in the case of ultra-short pulses. The analytical 
proof that those two solutions stay close, at least for an evolution up to 0(l/e) is 
a challenging task and subject to future research. 
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Figure 3. Comparison of the solution of Maxwell's equation and the 
short-pulse equation and the NLSE for an ultra-short pulse. Again, the 
figure above shows the initial pulse at a; = and the figure below com- 
pares the difference of nonlinear Maxwell's equation and the solution of 
the corresponding linear problem to the prediction of short-pulse equa- 
tion and NLSE. The parameters for this simulation are the same as in 
figure with the exception of 6 = 5 and = 20. 



Appendix A. Bounds on the roots of Q 

We begin by proving the results in Lemnias l3 . 3l and l3 .4l on the large p asymptotics 
of the roots of Q. We give all the details of the calculations needed to prove Lemma 
13. 31 and leave the very similar details of the proof of Lemma l3.4l to the reader. Begin 
by defining by 



4 = -|(r±v/?^^)+64 



(76) 
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Then Sj_ is a root of Q{s;p, e) if and only if ct^ is a solution of 

e'xo (-i(r± + 4 



2 



p2 
, 2 





p2 

If we let ^ = and define ^(cr^; ^, e) = g(cr° ;p, e) then we see that ^(0; 0, 0) = 
and 90-9(0; 0,0) = — ^JT"^ — Alu^ ^ 0, so the implicit function theorem implies that 
there exists a smooth function cr£(^,e) such that ^(^^(C, e); ^, e) = 0. Note that 
there exist constants A, B, C , and D which depend on F, f^nd xoi but not on e 
or p such that if we rearrange the equation g((T^(^, e); ^, e) = 0, (and replace ^ by 
it can be written as: 

|4(e, e)| < (1 - ^.Vp^)-^ + C|4(e, + e)t^ . (77) 

Since f7[|.(^ — 0, e = 0) = 0, (|77|l immediately implies that there exists Co,Ci > 
such that for p > Coc, 

I4(e,e)i<^, (78) 

which completes the proof of Lemma IX^ □ 
We now prove Lemma lS^ We first note that if p 7^ 0, Q has no purely imaginary 
eigenvalues. This follows by assuming that there exists such an eigenvalue - say 
s = ix, for a; G M. Inserting this into Q and equating real and imaginary parts we 
see that x must satisfy 

{eW,-x^)ip^-x')-e^x\o = (79) 
exr(p2 _ a;2) ^ . 

From the second of these equations we see that either a; = or a; = ±p. However, 
neither of these values of x solves the first equation and hence there are no pure 
imaginary roots. Next note that from Lemmas 13.31 and 13.41 we know that for p 
sufficiently large, all four roots lie in the left half plane. But since the roots vary 
continuously with p, the only way we could obtain a root with positive real part 
was is one of the roots passed through the imaginary axis. We have just seen that 
there are no pure imaginary roots for any non-zero value of p and hence we never 
have a root with positive real part. □ 
Finally, we prove Lemma 13.51 Note that if we define s — es and q — p/e, then 
Q{s;p,€) = Q{s,q,e), with Q{s,q,e) = + s^F + ujq)(s'^ + q^) + s^Xo, and the 
fact that the roots of Q can be written as s ^^{p) — es ^^{p/e) follows. Next note 
that for q small, an easy perturbative argument shows that Q has a pair of complex 
conjugate roots Sj. of size 0(1) (which correspond to the roots Sj_ of Q) and a pair 
of complex conjugate roots s\. — ztiq + O(g^). Thus, Lemma f3. 51 will follow if we 
can show that for all values of q the roots of Q are distinct. Since the coefficients of 
Q are real, the only way it can have multiple roots is if there is a multiple root on 
the (negative) real axis, or if there is a double complex root s and a second double 
root equal to the complex conjugate of s. We can immediately rule out the first 
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possibility by noting that 

Q{S, q, e) = s\s^ +Sr + uj^+ xo) + q^S^ + ST + oj^) . 

But since we have assumed that < 4wq, (s'^+sT+ujq) > and {s'^+sT+Uq+xo) > 
(for real values of s) and hence Q > for all real values of s. To rule out the 
possibility that the roots of Q are of the form s and s assume that that one can factor 
Q = {s — s)^(s — s)^. Expanding this expression and equating coefficients of like 
powers of s with the expression for Q above one finds a similar contradiction. □ 



Appendix B. Details on the derivation of NLSE by renormalization 

group method 

In this appendix we show how to solve for the higher orders approximating 
the solution of (|51() . Following the steps that led us to H58|) . let us now find the 
third order approximation. First, we need to collect 0{e^) terms. The usual way of 
obtaining these is simply plugging the previous ansatz u = euo+e'^ui+e^U2+- ■ ■ into 
(|51|l and collect appropriate terms. This will, however, lead to highly complicated 
algebraic calculation. We now approach this problem by making a different ansatz. 
Since we have already obtained the second order approximation of the solution, we 
assume that 

u ~ u^^^ + e'u2 + • • • . 

For the nonlinear part of the equation H51() we collect 0{e^) terms which will give 
rise to secular growth. To do this, we note that 



^U{X, tf = ^g3g3.(4.-^t) f ^ g ~ ^)3 

\ at J 

- e3e^(^^-"*) (^[b^ + 2i:bi^^ \U{x,t)\''U{x,t) 

_ g3g-»(/3.-it) |^3~2 ^ 2Cji^^ \U{x, t)\''U*{x, t) 

- e3g-3^(3x-i,t) J^g~2 ^ QCoi^^ U*{x,tf + 0{e^). (80) 



Thus, at the third order 0{e^), we have 
{£2 + /3')w2+^o(e±3^('^^-^*°)) 

+ (-/3/3"^A(a;i, ii) + 3ci2A(xi, ii)|A(a;i, e^(^---to) 

+ ^^A*ixi,h) + 3Cu^A*ixi,h)\Aixi,h)\^^ e-K/3.--to) = 0, 

where J^Qi^e^MPx-uito)-^ contains the terms proportional to Q^MP^-^tg) ^ ^^^^ A{x,ti) 
is from (|58|l . Note that the terms that will lead to resonances (secular growth) are 
the terms proportional to e^'^^^""^*). Solving the above differential equation, we 
find U2 = A2{ti)e"^^''~''^*°'> + P{x,to,ti) + complex conjugate, where ^2(^1) is a 
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function depending on the initial condition, and 

P{xMM) = -z:E(|^|^A(xi,ti)-^£A(a;i,ti)|A(xi,ti)ne*('^^-^*«) 



+ j£'o(e±^*(^^-'^*«)), 

where ^Q(^e^^^^fi^~'^^°'>) contains the terms proportional to Q^MPx-uito) ^ ^^^^ 
solve explicitly for the coefficients for these terms, and find these remain order 1 at 
all times. Now, the third order approximate solution is 



/3" 92 3 q2 



A2{ti)-ix — --2A(xi,ti) - -^A(xi,ti)|A(xi,ii)r e 



2 Mfix-Cjto) 



+ complex conjugate. 
Here, to get rid of the secular terms, we need to find Yj{x,ti) satisfying that 

S(a;, ti) \x=o = A(a;i,ti) +£2^2(^1), (81a) 

' = e'z(-^^S(a;,ii) + --^I](x,ti)|E(a;,ii)M . (81b) 

This leads us to introduce a new independent scale X2 — f^x in addition to xi — ex. 
Therefore, from H81|l . we finally obtain 

SU,=o - A{x,,ti)+e^A2{h), (82a) 
^ . ,Y4li^ + ^2!E|EFV (82b) 

dx2 y 2 dtf 2 /3 ' ' y ^ ' 

The equation Ij82b|) is the cubic nonlinear Schrddinger equation. We notice that the 
function E now depends on xi,X2,ti, i.e., E = I](xi, 0:2, ii). From here, together 
with (j57b|) . we obtain directly (|50|) by setting S(x, ti) = S(a;i, 2:2, ti). Following the 
similar steps, we can extend the results to the higher order, 0{e^) approximation. 
First, we assume the ansatz, u = u'-^^ + e^u^ + ■ ■ ■ and recall (|80|l . Recall dt — edti. 
Then, at the order 0{e^), we find 

(^ + Au, + e^(^--^*»)fi/3/3'''(-^)|^E + .f6^-3^^ 



\dx^ ' ' J yz"^'^ ^ 'dtl \ P dt 



e 



iircs 0; 

where .Fmcs contains the non-resonant terms, e.g. terms proportional to e='=3H/32:-wto) ^ 
g±4i(/3a;-wto) g^(,_ \hsX are generated when we insert U2 back into the nonlinearity. 
Solving the above differential equation, we find 

»3 = ^j(t,)e'<*--'">+jf„„. 



2i(3 

complex conjugate. 
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where ^3(^1) is a function depending on the initial condition, and J-'nrcs contains the 
terms proportional to higher powers of e^iP^-^^o) whose coefficients remain order 1 
regardless of x. Hence, we obtain the fourth order approximate solution, 

u<-^^ = eS(a;i, 2:2, tOe^t^-"-^*") + e3j£-o(e±3''(^--*o)) + e*^^,.,^ 

+ complex conjugate. 
In order to eliminate the secular term, we need to find V{x,ti) satisfying 

Vix,h)\,=o = S(2;i,a;2,ti) + e3v43(ti), (83a) 

This makes it natural to introduce a new scale 0:3 — e^x. From l)83b|l it follows that 



r = >'|i^— (6^-3^1 ^W). 
Since -^^^I^^P includes \ V\'^— — V and — we see that the Raman scattering 

uti ot\ Oti 

and the self-steepening terms appear at this high-order approximation. 
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